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Abstract 

We study the relaxational dynamics of the one-spin facihtated Ising model 

introduced by Predrickson and Andersen. We show the existence of a critical 

time which separates an initial regime in which the relaxation is exponentially 

fast and aging is absent from a regime in which relaxation becomes slow and 

aging effects are present. The presence of this fast exponential process and its 

associated critical time is in agreement with some recent experimental results 

on fragile glasses. 
02.70.Lq,75.10.Nr,64.60.Cn 
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I. INTRODUCTION 



The subject of glassy dynamics has received a lot of attention in the last years |l|]. During 
a fast enough cooling process real glasses do reach a metastable glassy phase of free energy 
higher than that of the crystal phase. Apparently the glass transition behaves as a purely 
kinetic phenomenon and the glass does not equilibrate when probed in a time scale smaller 
than the relaxation time. 

Laboratory experiments can measure one time extensive quantitites like enthalpy and 
its associated specific heat and also the two-times correlation function by measuring the 
scattering processes. These spectra give direct information about the relaxational processes 
which take place in glasses. One of the most studied relaxational processes in glasses is 
the so called structural or a relaxation which yields the structural relaxation time. While 
the a-relaxation is a slow process there are other faster processes which have been observed 
experimentally. Close to the glass transition two fast processes have been observed: 1) the 
/3-relaxation process predicted by the Mode Coupling Theory (MCT) [|| and observed in 
dielectric response measurements and 2) a faster process of order of picoseconds observed 
in neutron scattering experiments 0. In this last case, evidence has been reported on the 
existence of a critical time from the croossover from Debye (exponential in time and diffusive 
in space) to non-Debye relaxation |^. This critical time follows a temperature dependent 
Arrhenius behavior. The purpose of this work is to show that the existence of this critical 
time is an essential ingredient of some kinetic models with short range constrained dynamics. 

Several types of models have been proposed to understand the dynamical behavior of 
real glasses. All of them have in common the presence of a certain type of frustration. 
These models can be classified in two large classes. In the first class of models, there is 
frustration in the energy function. During its dynamical evolution these systems move in 
phase space avoiding configurations of higher free energy. The dynamics can be very slow 
due to the existence of energy barriers (strong metastability [||) or due to entropy barriers 
(strong marginality 0). Spin glasses [0 belong to these large class of systems where, in the 
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most general case, disorder is not essential and can be self-generated by the dynamics [§. 

While the first class of models (at least, in the mean-field approximation) seem to capture 
the experimentally observed features related to the (slow) a-relaxation process and the (fast) 



/3-relaxation process |Tl| it is still unclear how much they can account for this observed new 
type of fast Debye relaxational process [^]. 

In the second type of models the frustration is directly introduced in the dynamics. In 
this case the free energy landscape can be very simple but only certain transitions between 
configurations in phase space are allowed. These models are known under the name of 
Constrained Kinetic Models (CKM) a nice example being the n-spin facilitated Ising 
model (nSFM) introduced by Fredrickson and Andersen |[T0[| . 

There are few studies (theoretical as well as numerical) of this simple model but we 
think it contains some of the fundamental processes observed in real glasses. We will study 
the dynamical properties of one of the simplest models belonging to the aforementioned 
second class. In particular we will concentrate in the ISFM (to be defined below) at finite 
dimensions. We have observed that there exists a characteristic time t^, independent of the 
dimensionality of the system, below which relaxation is exponential and aging is absent and 
above which the relaxation becomes non-exponential and aging appears. This critical time 
follows an Arrhenius law with with the temperature and suggests a connection with some 
fast processes recently observed by the experimentalists 0. 

The paper is organized as follows. In the next section, we define the ISFM and the main 
observables we are interested in. Section three contains some exact relations for one-time 
staggered quantities at any dimension. This reveals the existence of two fast processes. 
Section four presents numerical simulations in the equilibrium regime and also in the off- 
equilibrium regime which evidence the existence of these fast processes in the two-times 
correlation and integrated response function. After the conclusions we present in the Ap- 
pendix the exact zero-temperature solution of the ISFM in one dimension. 
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II. THE ISFM MODEL: DEFINITION AND OBSERVABLES 



Let us take a set of field variables in a lattice of dimension D and a Hamiltonian 
H{(j){x)}. Let us consider an observable 0{t) which depends on time t through the configu- 
ration of the system O({0t(x)}). Now we define a discrete time dynamics for this system. In 
what follows we will consider a discrete Monte Carlo (MC) dynamics with random updating. 
A point of the lattice is randomly selected and a change of the variable (pix) is proposed. 
The rate variation of O in iV elementary moves (one Monte Carlo step (MCS)) is given by, 

= PiM^,t))WiM^,t) ^ (t>'tix,t))AO{t) (1) 

where AO{t) = 0{(f)[{x)) — 0{(f)t{x)) is the change in an elementary move of the set of 



fields (j)t{x), (■ ■ ■) stands for the average over all the possible transitions and P{(f)t{x)) is the 
probability of the configuration (f){x) at time t. 

We consider transition probabilities W of the form, 

W{Mx) <P[{x)) (X Mm(exp(-/3[if - 0^ (2) 

where the term a{(f)t, (p't) is temperature independent and cannot, in general, be absorbed 
in the energy function (note that the term a{(f)t, (p't) can be zero for certain transitions while 
this can never be the case at finite temperature for the MC dynamics). In the simplest case 
where a = 1 we recover the usual Metropolis algorithm for the Monte Carlo dynamics. Our 
purpose is to study a simple model where frustration only appears in the dynamics via the 
term a{(j)t,4't) ^'^ fact that some transitions between configurations are forbidden. 

A simple model of this type is given by the ra-spin facilitated Ising model (nSFM) in D 
dimensions |jlO[ . To each node of the lattice we attach a spin variable which can take the 



values 0, 1. The energy of the system is defined by the number of spins with value equal to 
one (with a minus sign), E = —Yl,i<^i- By defining the new set of variables s = 2(T — 1 we 
recognize in the previous expression the energy of an Ising paramagnet in a magnetic field 
h = ^. Fredrickson and Andersen proposed a constrained dynamics for the nSFM in the 
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following way: a randomly chosen spin in the lattice is selected and a flip of that spin is 
proposed (o"i —>■ a'^ = 1 — ai). This change is accepted if at least n of its nearest neighbours 
are and according to the probability Min{e~^'^^ , 1). The interest of this model is based on 
the fact that the non trivial (and glassy) dynamics is all contained in the term a{a{t), cr'{t)) 
(note that the Hamiltonian of this model has no interaction). Consequently, the dynamics of 
the nSFM is highly complex at inflnite temperature^ while in spin-glasses and other glassy 
models (where a{(f)t, ^^e dynamics in this limit is trivial. The nSFM model has been 

mainly studied in the case n = 1,2 |[10| , [T2|| . In this work we are mainly interested in the case 
n = 1 where some analytic results can be obtained. 



There are several physical interpretations of the nSFM [|T2|. The simplest one relates 
spin variables to the local compresibility of a fluid region. In this case regions with very 
high compressibility can facilitate the mobility of the neighboring ones while regions of low 
compressibility lead to a jamming of the dynamics. 

We deflne the set of variables Tj = 1 — a^. In terms of this set the transition probability 



W eq.m for the ISFM model reads P 



1 ^ 

W{a^ ^ a^) oc — exp(-/?a^) ^ (r^+e,, + r^-e J (3) 

^J.=l 

i.e. the transition probability at point x depends linearly on the magnetisation of the 
nearest neighbours. The set {e^; /i = 1,..,D} is a base for the the Z)-dimensional lattice. 
This transition probability satisfles detailed balance and is expected to generate an irre- 
ducible Markov process (in the termodynamic limit) in case n = 1 for any dimension at 
non- zero temperature (see the reference |12[ for a discussion on this point). 



While the thermodynamics of this model is trivial, its dynamics is much complicated 
and only partial results can be obtained in some cases, especially in one dimension. In this 
last case, the full dynamics can be exactly solved at zero temperature. Because we are not 
aware of this result in the literature we present it in the Appendix. Unfortunately we have 



^Note that at infinite temperature the system performs a non trivial random walk in phase space 
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not been able to close the dynamical equations at finite temperature. 

The general dynamical equation for an observable 0{t) in the ISFM model is given by 



dO{t) 1 



D 



AOix){-e~^ + (1 + e-^)r.)(^(r.+e, + r._ej) (4) 



dt 2D 

where AO{x) = 0{ax = 1) —0{crx = 0) stands for an elementary variation of the observable 
O at time t for the change of ax from to 1. Furthermore, the right-hand side has to be 
averaged over all points x of the lattice. 

We are interested in the energy-energy correlation function CE(t,t') (t' < t) defined as, 

where mT-(t) = Y^j^i Tiii) is the global magnetisation associated to the set of variables 
{Ti] i = 1, .., N}f^ Ce is normalized in such a way that CEit, t) = 1. 
We define the staggered one-point and two-point functions, 

^^o(t) = ^E(-l)^^-^"'V. (6) 
1=1 

CiW = ^ E E (-1)^^''-^''^ r.r.+e. (7) 

^^-^ u=l i=l 

where {i^; fi = 1, ..,D) are the coordinates of the point i in the D-dimensional lattice. The 
main interest of defining these quantities is that they can be exactly closed in any dimension 
(see next section). 

Another quantity of interest is the energy response function GE(t,t') defined by, 

t'<t (8) 

This function measures the change of energy of the system at time t if a small temperature 
change is done at t' < t. Eventually we will be also interested in the integrated response 
function, IRF(t',t) defined as 



^The same correlation function is obtained in terms of the variables {ai;i = 1, .., N} 
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IRF{t',t) = f'GE{t",t)dt" (9) 

J t 

In order to study the response of the system to a change of temperature it will be 
convenient to stay in the linear response regime. In this regime the change in energy is 
linear with the variation of f3. If the largest relaxation time (associated to the correlation 
function CE(t',t)) is strongly dependent with the temperature and if the perturbation of f3 
is not too small then we can expect strong deviations from the linear response regime. We 
will return to this point later. 

We should note that staggered functions can be defined only in finite dimensions and 
have no meaning in a mean-field version of the model. Our main interest is to show that 
some general results on the dynamics of the ISFM model can be inferred from the stag- 
gered one and two-point functions even though a complete solution of its dynamics at finite 
temperature is still lacking. 

III. FAST PROCESSES IN FINITE DIMENSIONS 

In the finite-dimensional case some general results can be derived for the staggered 
magnetisations eq.(|^). In the specific case of the ISFM the time evolution equations for 
Co(t), Ci{t) exactly close without the need of introducing hierarchies. Using eq.(^ it is easy 
to check that the time evolution equation for Co{t) reads, 

^ = -e-^C,it) (10) 

For Ci(t) the following result is obtained, 

^ = -(l + e-^)C,(t) (11) 

This simple result is a consequence of the particular dynamics for the ISFM as defined in 



eq.(^). In case of the 2SFM |]I2[ it is not clear if such type of relation exists. Unfortunately 
we have not been able to generalize this closure of equations in case of higher order corre- 
lation functions. It is important to note that staggered correlation functions are zero for 



homogeneous conditions (uniform or random for instance). Hence they cannot relax in time 
if initially they are at their equilibrium values. But it is reasonable to expect the response 
of the quantities Co and Ci to a staggered magnetic field to be also of exponential type. 

From previous equations (|1^) and (^TJ) we observe the existence of at least two charac- 
teristic time scales. One of these time scales is to = 1/(1 + e~^), the other one is tc = e^- 
The time scale to does not diverge at zero temperature and is nearly independent of the 
temperature (to = | at infinite temperature and to = 1 at zero temperature). This time 
scale corresponds to the relaxation time for the Ising paramagnet in absence of interaction. 
Hence, the time scale to corresponds to the relaxation of a single spin in the presence of a 
heat-bath and is independent of the nature of the constrained dynamics. 

The time scale tc is more interesting. It diverges at zero temperature and has an Arrhe- 
nius behavior. Moreover, this time scale corresponds to the smallest relaxation time of one 
spin interacting with its nearest-neighbors and is independent of the dimension of the lattice. 
Note that this diverging time scale is a direct consequence of the constrained dynamics in 
the system. Then, we conclude that tc is the first relevant time scale as a consequence of 
the cooperative phenomena which takes place in the ISFM . We will see later that this time 
scale is associated to a new fast relaxation process. 

In addition to these two time scales there is at least another one teq associated to the large 
time decay of the equilibrium two-times correlation function eq.(|^). This is the largest time 
scale and would correspond the a-relaxational process observed in glasses ^ Diagrammatic 



approximations done by Fredrickson and Andersen [jlO[ for the two-times correlation function 
eq.(^) show that this time scale coincides with tc. Our numerical investigation (in agreeement 
with previous work [T^) suggests that teq is much larger than tc and increases with P faster 



e^. But more detailed investigations are necessary in order to a better understanding of this 



^Wether in the ISFM there is also the P relaxational processes as predicted by the Mode CoupUng 
Theory M is still unclear 
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point. 

It is natural to think that the exponential relaxation processes described by equations 
(p!0|) and ([TT|) could be present also in the relaxation of non-staggered two times quantities 
like the energy-energy correlation function eq.(^) and the integrated response function eq.(y). 
Even though we do not have a rigorous proof of this assertion we have performed Monte 
Carlo numerical simulations for large lattices which clearly show that this is indeed the case. 
In what follows, we will focus our research in the physical consequences of the fast relaxation 
process described by tc- We will see the existence of a fast process of exponential character 
which manifests in equilibrium and off-equilibrium measurements. Unless otherwise stated, 
all results will refer to the ISFM case. 



IV. NUMERICAL RESULTS 

A. Numerical algorithm 

We have performed Monte Carlo numerical simulations of the dynamical equations in 
one and two dimensions. Starting from an initial configuration, the spins in the lattice are 
randomly selected and flipped according to the probability eq.(^). Simulations were done 
for relatively large lattice sizes = with L = 32000 in one dimension and L = 200 in 
two dimensions with periodic boundary conditions. In the range of times we are interested 
in (short time processes) finite-size corrections are certainly negligible and we have checked 
this is really the case. In order to test our previous exact results eq.(|lO]) and eq. fpUD , we show 
in the figure 1 the exponential decay of Co{t) at different temperatures as a function of the 
rescaled time t' = te~^ in one dimension starting from the periodic condition 11101110.... 



B. Equilibrium results 

In order to investigate the equilibrium properties, we should start from an initial fully 
thermalized configuration. There are several arguments which show that a dynamical transi- 
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tion is absent in the ISFM in finite dimensions [0. In this case, an equihbrium configuration 



can be easily built because the model has no interaction in the energy function, i.e. the cor- 
relation length is zero in equilibrium. We start from a random initial configuration with 
energy equal to the equilibrium energy, the system is let evolve and we compute the cor- 
relation function CE{t,t'). At equilibrium, the CE{t,t') is time-translational invariant and 
depends only on the difference of times CE(t, t') = Cj^it — t'). We have carefully checked this 
point. We have repeated this measurements at different temperatures. We have observed 
that correlation functions display two regimes separated by the critical time tc = . In the 
first initial regime (for times smaller than the critical time tc = e^) correlation functions 
decay exponentially fast. Empirically we find that the relaxation time associated to this 
exponential process is (with very high precision) 2tc, 

C|?(t)^ea;p(-t/(2g) = (Co(t))^ ; t « (12) 

In the second regime, for times larger than tc, the relaxation turns out to be much 
slower with stretched exponential behavior at high temperatures which becomes close to a 
power law behavior at low temperatures. We think that no conclusive relaxation behavior 
can be guessed numerically in this second regime. A crude estimate of the relaxation time 
teq for the second regime can be obtained (integrating the equilibrium correlation function 
respect to the time) and we obtain values larger than tc- This is an interesting point which 
deserves a more detailed investigation and should yield the relaxation time t^g for the slowest 
relaxational process^. Figure 2 shows C^{t) in one dimension at three diferent temperatures 
T = 0.15,0.2,0.40. The values of the critical time associated to these temperatures are 
tc — 780, 150, 12 respectively. These critical times are indicated with an arrow in the figure. 
The initial exponential decay eq.(|T2|) is also depicted in the figure as a continuous line. 



^Stretched exponential behavior already appears in the Glauber dynamics of the ferromagnetic 
one-dimensional Ising model. In this case it is difficult to guess the correct relaxational behavior 
from numerics. See the discussion of Brey and Prados on this point [Q. 
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Similar results have been obtained in two dimensions where we expect the this fast process 
to be independent on the dimension. This is shown in figure 3 where we plot the C|?(t) 
for two temperatures T = 0.15, 0.25 in one and two dimensions and the expected relaxation 
eq.(|T2|). The values of the critical time associated to these temperatures are tc — 780,55 
respectively. These times are indicated with an arrow in the figure. From this figure it can 
be clearly appreciated that the regime t < tc (the region of small times limited by the arrow) 
is independent of the dimension of the system as argued previously. The empirical relaxation 
eq.(l^) which fits pretty well the fast decay process is also depicted by a continuous line. In 
the second regime t > tc relaxation becomes dependent on the dimension and slower in one 
dimension than in two dimensions as expected (dynamical constraints in the ISFM tend to 
forbid less paths in phase space as the dimensionality of the lattice increases). 

We have also analyzed the integrated response function at equilibrium IRF{t',t) = 
IRF^'^{t — t'). In order to measure it we have prepared the system at equilibrium at a 
temperature (3. After some time we change the temperature of the system by the quantity 
A(3 and we let the system evolve at the new constant temperature f3 + A(3. If the change of 
temperature was applied at zero time then we have, 

/i?F-(t) = ^^--f^-^'^' (13) 

In order to be in the linear response regime it is necessary to make the change Aj3 small 
enough. But it is important to note that Af3 cannot be arbitrarily small, otherwise the 
response of the system to the temperature change is very small and the IRF'^'^{t) measure- 
ments become too much noisy. 

In figure 4 we plot the IRF^'^{t) normalized to its infinite time equilibrium value 
Ji?F^^(oo). This quantity converges to 1 in the inifinite time limit. Numerical experiments 
are shown in one dimension at two different temperatures T = 0.4 and T = 0.25 (data 
has been averaged over 20 different runs), for two different signs of the perturbation A/5. 
The results of the figure display some deviations from the linear response regime specially 
at large times and low temperatures. As commented previously, the equilibrium relaxation 
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time teq is strongly dependent on the temperature and increases with (3 faster than does tc- 
Then we expect the response of the system to a perturbation of the temperature to display 
departures from linearity in the part of the relaxation process controlled by the time scale 
teq. In the part of the relaxation process dominated by tc this effects should be smaller. 
Support in favour of this argument is shown in the figure 4 where the arrow indicates the 
value of the critical time for the two temperatures. The dependence on the sign of the per- 
turbation seems to be stronger for the slowest part of the relaxation (i.e. for times t > tc) 
than in the fast part of the relaxation (i.e. for times t < tc). In order to get more clear 
cut results it would be necessary to go to lower temperatures. Unfortunately the IRF^'^{t) 
becomes too much noisy even though we have been able to confirm this trend. 

C. Off-equilibrium results 

In this section we want to investigate the role of the critical time tc for non equilibrium 
relaxation processes. We are going to show that tc sets a minimum time scale above which 
non-equilibrium relaxations display aging effects. In the following we will use the notation in 
which t' ^ tw,t ty^ + t {tyj stands for waiting time). We have done numerical simulations 
of the ISFM in one and two dimensions measuring C Eitw^t^ + i). 

Results for the CE{tw,tw + t) are shown in figure 5 in the one dimensional case. We 
start from a random initial condition]^ with initial energy E{t = 0) = —0.5 quite far from 
its equilibrium value {E^''{[3) = — l/(l + e~^)). We computed the C£;(t^,t^ + t) for different 
values of the waiting time t^, = 10, 10^, 10^, lO'', 10^ at T = 0.1 (at this temperature tc — 
22000) in one and two dimensions (we show only the results in one dimension, in two 
dimensions they are qualitatively similar). In order to clearly appreciate the qualitative 
trend of the data we only show the values of t within the range 10'^ — 10^ MCS. Figure 



°The same results were obtained starting from a periodic initial configuration like for instance 
11101110.... where Co(t = 0) and Ci{t = 0) are different from zero. 
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5 shows two regimes. In the first regime {t^ = 10, 10^, 10^ < tc) aging effects are absent, 
i.e. the relaxation curve CEitwitw + t) only depends on t. In the regime = 10^ > tc 
aging effects appear (i.e. the full relaxation curve CEitw.tw + 1) also depends on tw). The 
continuous line is the equilibrium exponential behavior eg. (p!^) . 

As emerges from figure 5 the aging effects in the correlation funtion are very small even 
for tu, > tc at least in the region of times where the value of the correlation function is not 
too small. In fact, a scaling of the type CE{tw, tw+t) = f{t/tw) will not wotk at all. This is 
surprising since one would expect in the off-equilibrium slow regime a dynamical behavior 
plagued of strong aging effects in the correlation function. Possibly, this is a consequence 
of the particular correlation function used. In what follows we will see that the response 
function displays clear aging effects. 

We have measured the off-equilibrium integrated response function starting from a ran- 
dom initial configuration. The system is let evolve for a time at constant (inverse tem- 
perature) (3. At time t„, we make a copy of the system and we make it evolve at the constant 
new temperature f3 + A(3. After ty^ we measure the difference between the magnetisations 
of the two copies evolving with the same thermal noise but different temperatures. In this 
way we compute 

IRFit^, + t) = -^^^^^(' + ^-^--^^(^ + '-) (14) 

In figure 6 we show the normalized integrated response function IRF{tw,tw + 
t)//i?F^«(oo) at T = 0.1 (tc ^ 22000) for different values of = 10, 10^, 10^, 10^ 10^ 
in the one dimensional case. We can clearly appreciate the existence of the critical time tc 
which separates a regime where there is no t^i dependence from a regime where the system 
ages. As t^ increases the response of the system becomes smaller which is a typical feature of 
aging |T5|. Note that the equilibrium value of IRF'^'^{t) in the window of time shown in the 
figure 6 is practically zero (and converges to 1 for large times) while in the off-equilibrium 
regime the response reaches a (negative!) aproximate value of —10'^. Hence the response of 
the system is three orders of magnitude larger and negative in the off-equilibrium regime 
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than in the equihbrium case. 



V. CONCLUSIONS 

This work has been devoted to the study of one of the simplest constrained kinetic Ising 
models, i.e. the so called nSFM introduced by Friedrickson and Andersen in the special 
case of n = 1. This is an interesting model with an extremely simple energy function where 
the frustration is contained in the dynamical rules which forbid certain transitions between 
configurations in phase space. 

Our main interest has been the research of fast processes in the ISFM. In this model 
the rate of variation of the energy in one-point of the lattice is linearly coupled to the 
energy of the nearest neighbors. This is the simplest case one can consider, while the 2SFM 



corresponds to a quadratic coupling between the nearest neighbors sites in the lattice |]T0 
In case of the ISFM some exact closed dynamical equations can be obtained for the one 
and two-point staggered functions Cq and Ci in any dimensions. This reveals that both 
staggered functions Cq and Ci decay exponentially fast (with characteristic times tc and 
to) independently of the dimensionality of the lattice, to is the relaxation of a single spin 
uncoupled to its nearest neighbors embedded in a thermal bath (hence, independent of the 
dimension) and tc is an Arrhenius temperature dependent relaxation time resulting from 
the dynamical constraints. Then we expect the existence of these fast processes in other 
non-staggered quantities like the energy-energy correlation function and also the integrated 
response function. While we have not found a precise demonstration of this result we have 
given strong numerical support to this hypothesis by measuring the equilibrium and non- 
equilibrium behavior of the correlation function eq.(|^) and integrated response function 
eq.(^) in one and two dimensional lattices. In the equilibrium case the correlation function 
eq.(^) decays exponentially fast with an (empirically found) relaxation time 2tc in the regime 
of times t < tc while relaxation becomes slower in the regime t > tc. Furthermore, the 
integrated response function displays strong non-linear effects in the regime t > tc specially 
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at low temperatures. Concerning the off-equilibrium behavior we find that aging effects are 
absent for values of the waiting time less than the critical time. This is nicely observed in 
the behavior of the integrated response function (figure 6). 

Now we should discuss to what extent our results are general and not confined to the 
particular ISFM case. The equations (|TU|) and ( [TT| ) are exact results in the ISFM case and 
we expect they are not more valid in the nSFM with n larger than 1. We have performed 
some numerical simulations in case of the 2SFM but we have not found evidence on the 
existence of this critical time. This could explain why some previous numerical works on the 



2SFM [T^ were able to fit the equilibrium relaxation functions to a strecthed exponential 
behavior while in the ISFM this is not possible due to the existence of two different regimes. 
In numerical simulations one is able to explore only relatively small scales of time. It is clear 
that an initial exponential process would make numerical fits of the slow regime difficult, 
specially in the low temperature phase where the critical time starts to become large. 

Also in the realm of disordered systems (for instance in spin glasses ||^) we are not aware 
of the existence of this critical time, possibly because in that case frustration is directly 
introduced in the energy function and not in the dynamics. We are tempted to conclude that 
(1) the existence of a temperature activated critical time and (2) the presence of exponential 
decay processes in the relaxation of some correlation functions in the short-time regime, are 
both strictly related to the nature of this type of constrained short-range dynamics. 

Summarizing, we have found evidence about the existence of two fast processes in the 
ISFM model. The first process with characteristic time to corresponds to the relaxation of 
single spin in the lattice. This is a trivial process not related to any cooperative effect in 
the lattice. The second process is a fast exponential process consequence of the dynamical 
constraints. Physically it would correspond to the relaxation of some coupled degrees of 
freedom of the system. The results which emerge from the study of the ISFM are in 
agreement with some recent experimental findings by Colmenero et al. on fragile polymer 
glasses. By doing neutron scattering measurements in glasses they claim on the existence of 
a temperature activated critical time (with small energy barrier) which separates two well 
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defined time regions. In the first region the relaxation is exponential while in the second 
regime relaxation is much slower with higher relaxation time. This feature seems to be 
captured by the present model. 

It would be very interesting to analitically solve the dynamics at finite temperature of 
this model (at least in one dimension) in order to confirm the numerical findings of this 
work. This would also shed light on the existence of other relaxational processes (like the a 
and /? relaxation) as predicted by the Mode Coupling Theory. 
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VI. APPENDIX 

In this appendix we present the exact solution of the ISFM in one dimension at zero 
temperature. Lets take a chain of spins ai and we define the new set of variables = X—oi. 
In terms of this set of variables we define the following set of correlation functions, 

1 ^ 

^A.(t) = T7 E 7-r(t)r.+i(t) ... r,+fc(t) (15) 

The magnetisation m,- = j^Yl!^=\Tr = -Do(t) is the first term of this hierarchy. Using 
eq.(§) we can derive the time evolution of the Dk(t) at zero temperature. We get, 

^ = -Dk+i-kDk] k>Q (16) 

We introduce the generating function 

G{x,t) = j:T[Mt) (17) 

fc=l 

In terms of this generating function we have the following partial differential equation. 
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^^ = -(l + x)^^ (18) 

This linear partial differential equations is readily solved yielding, 

G{x,t) = Goiil + x)e-'-l) (19) 

where Go{x) = G{x,0) is the initial condition. We can obtain the different set of mo- 
ments, 

In particular we get for the magnetisation m,- = Go{exp{—t) — 1). In the large time 
limit it converges to Gq{—1) depending on the initial condition. For the particular initial 
condition cij = 0, Go{x) = and the magnetisation rrir does not converge to its equilibrium 
value (m^"^ = 1) but to 1/e. 

We can also compute the two-times correlation function eq.(^). At zero temperature we 
find, 

_ ^ - ^(^) _ ^rjt) _ Goje-* - 1) 
^^^^ ' - 1 - m(f ) - - Go(e-*' - 1) ^^^^ 

In figure 7 we compare previous equation ( pi] ) with the numerical results for different 
values of t'. 
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Figure Captions 

Fig. 1 Exponential relaxation of Co{t) in one dimension for several temperatures T = 
0.1,0.2,0.3,0.4,0.5. 

Fig. 2 C^lt) in one dimension at temperatures T = 0.15,0.25,0.4. The continuous lines are 
the exponential relaxations eq.(|T2p. 

Fig. 3 Cj^it) at temperatures T = 0.15 {D = 2squares, D = lacrosses)) T = 0.25 {D = 
2{times),D = l{circles)) in one and two dimensions. The continuous lines are the 
exponential relaxations eq.(|I^). 

Fig. 4 Normalized IRF'^'^it) in one dimension at temperatures T = 0.4 with A(3 = 
0.2{squares), —0.2{crosses) and T = 0.25 with A/3 = OA^rhombs), —OA^dots) 

Fig. 5 CE{ty,, ty, + 1) for different waiting times = 10, 10^, 10^ 10^, 10^ in one dimension at 
T = 0.1. For tyj = 10, 10^, 10^ the relaxation curves superimpose and aging is absent. 

Fig. 6 Normalized IRF{t^,t^ + t) at T = 0.1 for several values of and A(3 = 0.5. The 
symbols are guide to the eyes. Note that IRF^'^{t) is nearly zero in this range of times. 

Fig. 7 CEitw,tw +t) for different waiting times tw = 0.0125,0.125,1.25 at zero temperature 
in one dimension. The lines are the exact solution (21). 
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